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ABSTRACT 


Two mathematical representations of noise due to atmospheric turbulence are 
presented. These representations are derived and used in computer simulations of 
the Bartlett Estimate implementation of beamforming. Beamforming is an array 
processing technique employing an array of acoustic sensors used to determine the 
bearing of an acoustic source. Atmospheric wind conditions introduce noise into 
the beamformer output. Consequently, the accuracy of the process is degraded and 
the bearing of the acoustic source is falsely indicated or impossible to determine. 
The two representations of noise presented here are intended to quantify the effects 
of mean wind passing over the array of sensors and to correct for these effects. The 
first noise model is an idealized case. The effect of the mean wind is incorporated as 
a change in the propagation velocity of the acoustic wave. This yields an effective 
phase shift applied to each term of the the spatial correlation matrix in the Bartlett 
Estimate. The resultant error caused by this model can be corrected in closed form 
in the beamforming algorithm. The second noise model acts to change the true 
'^direction^of propagation at the beginning of the beamforming process. A closed 
form correction for this model is not available. Efforts to derive effective means 
to reduce the contributions of the noise have not been successful. In either case, 



Ill 


the maximum error introduced by the wind is a beam shift of approximately three 
degrees. That is, the bearing of the acoustic source is indicated at a point a few 
degrees from the true bearing location. These effects are not quite as pronounced 
as those seen in experimental results. Sidelobes are false indications of acoustic 
sources in the beamformer output away from the true bearing angle. The sidelobes 
that are observed in experimental results are not caused by these noise models. 
The effects of mean wind passing over the sensor array as modeled here do not 
alter the beamformer output as significantly as expected. More research is required 
to determine the cause of the sidelobe errors. Though smaller than expected, the 
effects of mean wind are quantified. If the wind may be represented by the first 
model, the effects of the wind may be eliminated. Correction for error caused by 
the second model is not currently possible. 
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CHAPTER I 
INTRODUCTION 

Beamforming is a signal processing technique which uses an array of acoustic 
sensors to determine the bearing of an acoustic source. The outputs of the sensors 
are summed with weights and time delays to form a beam. The time delays are 
called steering vectors and their function is to sweep the beam through a range 
of assumed bearing angles. When the assumed direction corresponds to the true 
bearing angle, the magnitude of the beam is maximized. Beamforming is a passive 
detection technique. Data are digitally sampled from the outputs of acoustic sensors 
and then processed by computer. No energy is emitted as in active techniques such 

as radar. 

Technology employing the acoustic principles on which beamforming is based 
was used as far back as World War I. These principles will be discussed in detail 
in Chapter II. Such processes were first developed by the British and the French 
to detect approaching bombers and zeppelins in the air and submarines in the 
ocean. However, these methods did not have sufficient accuracy to localize targets 
for aiming weapons. Towards the end of World War II, technology was developed 
to use radar for determining the range and bearing of enemy craft. Since World 
War II, most uses of passive detection techniques have been employed in water 
for the purpose of tracking. There are also geophysical applications for gathering 
information on seismic activity as discussed by Kleyn (1983). These applications 
are well documented and much research has been done in these areas. 
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Recently, interest has been renewed in the application of beamforming tech- 
niques to the tracking of acoustic sources (particularly helicopters) in the atmo 
sphere. The advantage of beamforming is its passive nature. The acoustic sensor 
array emits no energy aimed at the target that the target can detect. This is not 
the case with the use of radar which emits a beam of electromagnetic energy and 

records echoes as the beam bounces off potential targets. 

The accuracy of the beamforming process is degraded by noise in the signal 
from the acoustic sensors. The major portion of the noise is introduced as 
the acoustic waves pass through the medium. Variations in temperature and 
velocity of the air or water cause the path and speed of the acoustic waves to 
be varied. The result is increased uncertainty of the bearing angle of the target 
source. False indications may also be present. Due to the widespread use of 
beamforming applications in the ocean, the problem of noise in water has been 
addressed considerably and its characteristics are fairly well understood as discussed 
by Burdic (1984). The nature of noise due to atmospheric turbulence and how noise 
affects the beamforming process for the purpose of tracking in air has not received 

much attention. These topics are the subject of this report. 

The Bartlett Estimate is the beamforming implementation discussed in this 
report. The Bartlett Estimate is but one beamforming algorithm. There are other 
high resolution techniques such as the Maximum Likelihood Method first reported 
by Capon (1969) and the Linear Predictive Estimate presented by Johnson (1982). 
These high resolution algorithms outperform the Bartlett Estimate in some cases. 
For example, these techniques are much better at distinguishing between two 
separate sources at close bearing angles where the Bartlett Estimate might identify 
the two as only being one source. A major drawback of such high resolution methods 
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is that a matrix inversion is required. In cases like tracking where averaging time 
is short, the matrices to be inverted may be singular and the methods fail. The 
Bartlett Estimate does not require matrix inversion and therefore can be applied to 

tracking. 

Beamforming methods typically assume that noise is ideal. That is, the noise is 
assumed to be uncorrelated from sensor to sensor, spatially white and homogeneous. 
Therefore, the effects of the noise can theoretically be eliminated by long time 
averaging of the signal. Long time averaging is possible in cases where the source 
motion is small over the averaging time. However, in the case of a rapidly moving 
source or any case in which the averaging time is limited (e.g. tracking), the effects 
of noise are not negligible and cannot be averaged out. If a suitable representation 
of the noise due to atmospheric turbulence could be developed, it could then be 
used to modify the original data containing the acoustic signal and noise to yield 
a “cleaner” signal with which to carry out the beamforming process. Therefore, 
the accuracy of the beamformer would be enhanced. This report presents two 
such representations (models) for noise due to atmospheric turbulence. Computer 
simulations are performed to quantify the effects of noise having the form of these 
two models on the accuracy of the Bartlett Estimate beamformer. Both models 
require knowledge of the local wind conditions (speed and direction). The first 
model is idealized and a closed form solution is allowed. That is, the error induced 
by the wind can be eliminated within the beamforming algorithm. The second 
model is less idealized and a closed form solution is not available. Efforts to correct 


this error were unsuccessful. 
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CHAPTER II 

THE THEORY OF BEAMFORMING 

The primary objective of beamforming is to determine the bearing of an 
acoustic source. An array of acoustic sensors is used to gather data on the actoustic 
energy in the air. The sensors record the acoustic vibrations in the air. These 
vibrations include contributions coming from the source at some unknown bearing 
angle and contributions from noise present in the atmosphere. The output of each 
sensor is multiplied by a weight a, and given a time delay Wi. The adjusted outputs 
are then added to form the beam. Two assumptions are made: 

1) The sound waves passing over the array are plane waves. This assumption 
provides that the source is sufficiently far away that the curvature of the acoustic 
wave is negligible over the sensor array. 

2) The acoustic signal is only a function of time. No attenuation occurs over the 
distance traveled between sensors. The signal that each sensor receives is identical 
except for a phase shift (essentially a time delay) that arises from the different 
distances that the wave must travel from sensor to sensor. 

- Beamforming attempts to offset the phase shift (time delay) exactly at each 
sensor. When this is achieved, all the weighted sensor outputs are in phase and 
the magnitude of the beam is maximized. A beamformer is an example of a spatial 
filter. That is, acoustic energy arriving at bearing angles other than the true bearing 
angle are attenuated. The Bartlett Estimate studied in this thesis is presented by 
Nielsen (1989). The derivation of the Bartlett Estimate is presented next. 

A target acoustic source is located at some unknown bearing angle 9. A unit 
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vector in the direction of the acoustic waves propagating from the acoustic source is 
u. An array of M acoustic sensors is assembled and a coordinate system is defined. 
Each sensor has position vector f, and a reference sensor is designated. See Figure 1. 
Note that the unit vector -u is located by the bearing angle 9. Let the reference 
sensor be sensor 1. The speed of sound in the air is c. The output at the ith sensor 

is given by 

*,(<) = ,(< + + »<(<) (2-1) 

where s is the acoustic signal from the source and n is noise in the sensor output. 
The differential travel distance between sensor 1 and sensor i is r, • u and is shown 

in Figure 2. 

The time delay t, due to the differential travel distance is given by 

= ! l 'JL ( 2 . 2 ) 


This is the time the plane wave takes to travel between sensor 1 to sensor i. 

A bearing angle 9 0 is assumed and the unit vector pointing from this bearing 
direction is uq. The offsetting time delays r t0 are then calculated as 


t;o = ^ ( 2 - 3 ) 

c 

The assumed bearing angle is swept through a range of assumed bearing directions. 

The sensor outputs are given the offsetting time delays, multiplied by weights 
a t and summed. This sum is called the “beam” or “beamformer.” The beamformer 
output is denoted by g(t ) and is given by 


g(t) = 


M 

E 

i = 1 


0 L{ s(t 


+ 


u 0 


u 0 




(2.4) 
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When the assumed bearing angle corresponds to the true bearing angle, the 
offsetting time delays exactly cancel the delays of Equation 2.2, all of the sensor 
outputs are in phase and their sum is maximized. Assuming unity weights 
(a, = 1.0) and letting r i0 denote the offsetting time delays at each sensor as 

in Equation 2.3, the beamformer output is given as 
M 

g(t) = Ms(t) + £ n(t - Tj 0 ) (2.5) 

i = l 

for the correct value of offsetting time delay r; 0 . 

A. Frequency Domain Analysis 

Further insight into the function of beamforming is gained by transforming 
these expressions into the frequency domain. The application of digital computers 
to beamforming in the frequency domain allows the use of fast algorithms such as 
the Fast Fourier Transform (FFT). Taking the Fourier Transform of Equation 2.4 
yields the transform of the beamformer G(f , uq). G(f , uo) is a function of frequency, 
/, and assumed direction, uo, and is given by 

G(/,&.) = (2.6) 

i = 1 ^ ' 

where j is y/^1 and Xi(f) is the Fourier Transform of the sensor output 

X l (f) = S(f) exp | ~~ | + Ni(f) (2.7) 

Note that the offsetting time delays in Equation 2.4 are now represented by 
multiplication of complex exponentials which yield phase shifts in the frequency 


domain. 


10 


G(f,u 0 ) can be represented by a dot product of two vectors w and z which 
have components 

Ci = exp( — — — 1— — — I and i, = X,(f) 

The terms w t are the steering vectors. Note that they include the assumed bearing 
angle 0 O . Then the beamformer is 

G{f,u o) = w + -x (2-8) 

where + denotes the complex conjugate transpose. 

The power of the beam is a function of the square of the magnitude of the 
beam. With assumed direction u 0 and at frequency / the power is given by 

P(/,u o) = E{\G(f,uo)\ 2 } ( 2 - 9fl ) 

= w* E{xx + (2.96) 

= w + Rw (2.9c) 

where £{•} is the expectation ‘operator which yields the average value of the 
argument. R is the spatial correlation matrix. The fcnth element of the spatial 
correlation matrix is given by 

R kn = E{[S(f) exp{j27r/ } + iV fc (/)]x 

[S’(/)exp{-j27r/^-^} + Kif)}) ( 2 ‘ 10 ) 

where * denotes the complex conjugate. It is assumed that the signal and 
noise portions of the sensor outputs are completely uncorrelated (statistically 
independent). Therefore, the average of cross terms involving noise terms multiplied 
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by signal terms in Equation 2.10 are equal to zero. The spatial correlation matrix 
can then be expressed as 

*. = e{s(/)S-(/)« cp{;2,/ (Ft -^ }} - E{.WM)} (2.11) 

Defining the spatial correlation matrix for noise as 

_ E{n k {f)n m n {f)} (2A 2) 

Q^n = ~2 v 


the beam output power is 


P(/,u 0 ) = exp{j'27r/ 


(r fc - r-n) • u 




crlw' 


Qw 


(2.13) 


The expression in Equation 2.13 is the Bartlett Estimate for the power of the 
beamformer output. o 2 and are the intensities of the signal portion and the 
noise portion of the sensor outputs respectively. a\ is equal to £{5(/)S (/)}• 
The signal-to- noise ratio, SNR , is given by <r 2 J<r 2 n . Note that Q is normalized by 
a -2 . The beamformer scans the range of assumed bearing angles for each frequency 
component in the range of frequencies of interest. In practice, much computer 
processing capability is required for the purpose of tracking since these ranges must 
all be covered in a time period that is approximately one second. 

In Equation 2.13, a], <r 2 n , u and Q are unknown quantities. The SNR in 
practice is usually small (on the order of 0 dB). The relative strengths of the signal 
and noise are important in determining the ability of a particular beamforming 
algorithm to pick out the signal when it is embedded noise. Q is the chief source 
of uncertainty because its composition is unknown. Many idealized models have 
been proposed to try to quantify the nature of Q. The short averaging time 
characteristic to tracking causes poor performance of these models when applied to 
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the beamforming process. One such model is briefly discussed in the next chapter 
and three others are included in an appendix. 

In experimental results, three primary effects on the beamformer performance 
are observed to be caused by atmospheric conditions. These effects are as follows: 

1) Beam Shift - The indicated source bearing angle is displaced from the true bearing 
angle. 

2) Beam Widening - The main peak in the beamformer output at the true bearing 
angle is widened. This adds uncertainty to the actual location of the maxima in 
beam power in practical implementation. The beam width is measured at the point 
where the beamformer power magnitude is 50% of the maximum beam power. 

3) Presence of Sidelobes - Sidelobes are maxima in the beamformer power output at 
assumed bearing angles away from the true bearing angle. The sidelobes indicate 
that there are acoustic sources at other bearing angles when in fact there are none. 

The matrix Q from Equation 2.13 contains the contribution of noise due 
to turbulence in the atmosphere. This quantity is unknown. If a suitable 
representation of Q could be determined, it could be subtracted from the right 
hand side of Equation 2.13 and a more accurate beamformer output would result. 
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CHAPTER III 

APPLICATION OF BARTLETT ESTIMATE 
BEAMFORMER TO CURRENT WORK 

A. Array Geometries and Spatial Aliasing 

The range of frequencies over which a given geometry may be effectively used 
in beamforming is limited. The upper end of this range is limited by spatial 
aliasing. For a linear array with equally spaced sensors, the critical upper limit 
is the frequency whose wavelength A is twice the sensor spacing d. When the 
frequency is greater than this critical value, the actual direction of propagation u 
may be confused with other assumed values of propagation ilq. Aliasing causes 
considerable sidelobes thereby indicating sources where there are in fact none. The 
location(s) where spatial aliasing will cause false indications is uncertain. This is 
related to array geometry, the sensor spacing, d, the frequency, /, and the true 
bearing angle, 6. In practical applications, spatial aliasing is eliminated by low-pass 




Figure 3. - Linear Sensor Array 
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filtering. It is important to know the critical frequency for the given array geometry. 
For the linear array with a sensor spacing of 1.8288 m (6 ft), the upper frequency 
limit is 94 Hz. Figure 4 gives an example of spatial aliasing. Beamformer outputs 
for M = 10, d = 1.8288 m (6 ft), 9 = 45° and pure signal (no noise) are given for 
frequencies of 90 Hz and 140 Hz. The output for 140 Hz indicates a false source at 
129°. Spatial aliasing is a phenomenon characteristic to beamforming and occurs 
given these conditions regardless of the noise model used or the array geometry. 
However, it can be minimized by array geometry. For example, inner sensors of the 
nested triangular array may be used to process higher frequencies than the outside 
sensors. 

The lower range of frequencies is limited in that, as the frequency decreases, 

A increases and becomes larger in comparison to the sensor spacing. The result is 
a broadening of the peak in beamformer output at the true angle of propagation 
and therefore the true angle of propagation becomes more uncertain. Figure 5 gives 
examples of this characteristic. Beamformer simulations for 30, 45, and 90 Hz with 
pure signal output for equal values of d, M and 6 are shown. At 90 Hz, the half 
power point beam width is 15 degrees. At 45 Hz, it is 31 degrees. At 30 Hz, the 
half power point beam width is 54 degrees. The lower frequency limit should be 
set depending on the permissible beamwidth for a particular application. An upper 
limit recommended by Gerhold (1990) is the frequency whose wavelength A = Ad. 

Another consideration concerning the linear array configuration is symmetry. 
The output of a beamformer using a linear array is symmetric about the axis of 
the array line. That is, if the line connecting the sensors is the reference for zero 
degrees, any source, at 10° for example, will have an identical and yet false indication 
at -10°. This characteristic of the linear array is illustrated in Figure 6. In the 
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Figure 4. - Bartlett Estimate beamformer power output showing example of spatial 
aliasing. False indication is shown at 129 . 


EFFECT OF FREQUENCY ON BEAMWIDTH 
pure signal, d = 1.8288 m, M = 10, 6 = 45 



Figure 5. - Bartlett Estimate beamformer power output showing the effect of 

frequency, /, on the half power point beamwidth. 
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Figure 6. - Illustration of symmetry for the linear array of equally spaced acoustic 
sensors. True source is at 45°. False indication shown at 45 . 
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interest of space and clarity, most plots of beamformer output in this report will 
have 0° to 180° as a range for the assumed bearing angle 0 O - 


B. Description of Beamformer Simulations 

A linear array with equally spaced sensors as in Figure 3 was modeled. The 
sensor spacing d was 1.8288 m (6 ft). The half power point beam width is related 

to the number of sensors M. This is illustrated in Figure 7. 

For identical values of d, f and .6, beamformer simulations were performed for 
values of M equal to 10, 15 and 20. For M = 10, the half power point beam width 
is 15 degrees. For M = 15, it is 10 degrees. For M = 20 it is 9 degrees. Note that 
as M is increased, the beam power increases as well. The number of sensors used 
here in the main simulations was M = 10. The speed of sound, c, used was 344 m/s 
(1128.0 ft/sec). This speed corresponds to a temperature of 294° K (70° F) at sea 
level. Discretized sine waves at a frequencies of 30, 45 and 90 Hz were generated and 
taken as the signal input to the sensors. (Recall that the frequency domain analysis 
of beamforming processes one frequency component at a time.) These frequency 
values yield results that are representative of beamformer output over the useful 



EFFECT OF M ON BANDWIDTH 
d = 1.8288 m, 5 = 45°, f= 90 Hz 



Figure 7. - Illustration of the effect of the number of sensors, M , on the half power 
point beamwidth of the Bartlett Estimate power output. 
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range of frequencies for the 1.8288 m (6 ft) sensor spacing. This useful range is 
the range of frequencies that the beamformer can be used to process. This useful 
range is limited on the high end by spatial aliasing. The low end is limited by beam 
widening. This range is related to the ratio of d/X. Since the array geometries are 
typically fixed, the frequencies were varied while d was held constant to demonstrate 
beamformer behavior as the ratio d/X changed. The upper aliasing limit where the 

wavelength d = A/2 for d = 1.8288 m (6 ft) is 94 Hz. 

It was assumed that the electronics introduce no noise into the sensor outputs. 
That is, the sensor outputs were modeled as pure, discrete sinusoids. A 512-pomt 
Fast Fourier Transform was taken of the signal outputs to transfer the data into the 

frequency domain. 

All plots included to illustrate results present the magnitude of the beamformer 
power output on a linear scale. Typically, a decibel scale is used. However, a linear 
scale better shows the characteristics that are to be pointed out as results in this 

thesis. 

Simulations were performed for different values of 9 over the range of 0 to 
180° for 9 0 . Increments of one degree were taken. The frequency was 90 Hz. The 
beam width is dependent on the value of 6. As 9 approaches the extremes of this 

range, the beam becomes wider. See Figure 8. 

This beam widening is symmetric about 90”. For 0 = 0”, the half power point 
beamwidth is 50 degrees. For 9 = 30° it is 23 degrees. For 9 = 90°, it is 10 
degrees. Recall that the beamformer output is symmetric about the line connecting 
the sensors in the linear array. This line is at zero degrees. As the beam moves 
toward zero degrees from the right, there is a corresponding beam moving toward 
zero from the left and these two begin to join together as they approach zero degrees 



Figure 8. - Illustration of the effect of the true bearing angle, 0 , on the half power 
point beamwidth of the Bartlett Estimate power output. 
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(or 180 degrees). This contributes to a very large half power point beam width n 
the extremes. 

C. Criteria Used to Judge Performance of Noise Models 

The power output of the Bartlett Estimate beamformer and the effects of the 
two proposed noise models on it are simulated in this report. Three criteria are 
used in judging the performance of the noise models in comparison to beamformer 
simulations with no noise included. 

1 ) Shift in the location of the maximum beam power. This is the error in dete g 

the “exact” true bearing angle of the acoustic source. See Figure 9 

illustration of this shift. 

2) Half power point beam width. This is the width m degrees of the peak in 
beamformer output around the indicated true bearing angle corresponding to a 
power value of 50% of the peak power value. See Figure 10 for an illustration of 

the half power point beam width. 

3) Presence of sidelobes. The presence of sidelobes indicate that a particular model 
causes false indications. That is, that the model would show that there are acoustic 
sources at bearing angles where there are in fact none. See Figure 11 tor an example 

of sidelobes. 


ILLUSTRATION OF BEAM SHIFT 
pure signal, d = 1.8288 m, M = 10, f = 


^ beam shift m 
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Figure 9. - Illustration of beam shift in the Bartlett Estimate beamformer power 
output. A beam shift of 5° is shown. . 
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ILLUSTRATION OF HALF POWER POINT BEAM WIDTH 
pure signal, d = 1.8288 m, M = 10, 8 = 45 , f — 90 Hz 
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Figure 10. - Illustration of the beam width at the half power point of the Bartlett 
Estimate beamformer. 










Figure 11- - Example of the presence of sidelobes in the Bartlett Estimate 

beamformer power output. The true bearing angle 9 is 67°. A false indication 
of an acoustic source is present at 113° due to the sidelobe at that bearing angle. 
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CHAPTER IV 

TWO REPRESENTATIONS OF NOISE DUE TO 
ATMOSPHERIC TURBULENCE 

This chapter will present the two actual noise models proposed in the current 
work. The first is based on alteration of the speed of propagation and the second is 
based on shifting the direction of the propagation vector. The first is idealized and 
is correctable. The second is less idealized and efforts to correct it were unsuccessful. 

A. Modification of Bartlett Estimate 

A departure is taken from the Bartlett Estimate expression of Equation 2.13. 
The spatial correlation matrix for noise, Q , is omitted and beam power is calculated 
from pure signal 

P{f,u o) = w+Rsw 

Rs is an adjusted spatial correlation matrix containing only signal information and 
no noise information. This is obtained from Equation 2.11. The knth element of 
Rs is given by 

= S(f)S'(f) exp{j2ir/ (rfe ~ c rn) — } ( 4 - 2 ) 

For the linear array, the dot product in Equation 4.2 is given by 

( ffc - F n ) • u = ( k - n)dcos6 (4-3) 

Rather than account for the presence of noise by using the matrix Q, the noise 
is represented by an alteration of the phase shift information of Rs> H° w ^is is 
accomplished for each of the two noise models will be discussed in detail later in 
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this chapter. A wind speed of 13.4 m/s (44 ft/sec) was used. On flat land, this 
is a realistic order of magnitude mean wind velocity for a short sustained gust as 
presented by Stull (1988) and recommended by North (1990). 

Each case of modeling the noise due to atmospheric turbulence is compared 
to the calculation of beam power with no noise as given by Equation 4.1. This 
allows the full effect of the wind to be isolated and observed. 


B. Noise Model 1 


The first noise model assumes that the wind alters the speed of propagation 
of the acoustic waves. The direction of the propagation vector is assumed to be 
unaffected. Doppler effects are also ignored. The projection of the wind velocity 
onto the vector of the true propagation velocity adds to the speed of sound c 
and effectively changes the speed at which the plane wave travels through the 
atmosphere. The change in velocity alters the time delay which occurs as the 
wave passes from one sensor to another. Recall from Equation 2.2 that the time 
delay is dependent on the propagation velocity. Since the beamformer attempts to 
negate these phase shifts based on the propagation speed in still air, an error will 
be introduced. This error will result when the phase shifts represented in Rs and 
in the steering vectors w x will be equivalent but the assumed angle of propagation 

and the true angle of propagation will not be the same. 

Recall the definition of the fcnth element of the spatial correlation matrix for 

pure signal Rs 


Rs.fcn 


5(/)5-(/)exp{;27r/ 


(ffc ~ r n) 


C 


(4.4) 
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A unit vector pointing from the true bearing direction is u. The true bearing of 
the acoustic source is 0. Now consider a wind blowing at constant velocity V m the 
direction defined by a unit vector p. The bearing angle of the direction that the 
wind is coming from is 0 W . See Figure 12. Unit vector -p points in this direction. 
The wind acts to change the phase shift between two given sensors. 

Consider the time delays between the *th and nth sensors. The extra distance 
that the plane wave must travel after passing sensor n before sensor k receives it is 

Pkn given by 

Pkn = (n - k)dcosd ( 4 - 5 ) 

This distance is illustrated in Figure 13. The speed of sound in still air is c and 
then the corresponding time delay Tfc n is pkn/c- 

(n — k)d cos# (4.6) 

"Hen = ~ 

c 

Now let the wind velocity component along u be added to c. The adjusted 
speed of propagation is given by 

c -r V cos(#w — 9) 


Figure 14 illustrates the vector operations that yield this result. 

- The adjusted time delay which occurs now between the nth and fcth sensors is 


r a , fcn and is given by 


(n — k)d cos 0 

Ta ' kn = c + V cos(^w - O) 


The difference between time delays without wind and with wind (the difference 
between Equation 4.6 and Equation 4.7) is Td, n k and is given by 


Td t kn =(n - k)d cos 0 { - - 


c + V cos(Qw - 0) 
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(d) 

Figure 14. - Vector additions of wind and sound propagation vectors for Noise 
Model 1. (a) Wind vector, (b) Sound propagation vector in still air. (c) Projection 
of the wind vector in the direction of propagation vector shown in (b). (d) Resultant 
adjusted propagation vector with unaffected direction but with adjusted magnitude. 
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In the frequency domain, this time delay corresponds to a phase shift <t> which, in 
radians, is 2ir fr d M- A phase shift is represented by a complex multiplication in the 
frequency domain. This phase shift can be used to alter the phase characteristics 
of the matrix Rs by multiplying the knth element of Rs by e . This defines the 
spatial correlation matrix for the first model, R\. This phase shift is given by 

**» = 2 *H n - *)<'“’ ” c 7 v cos(0w - T ) ) (4 ' 9) 

The knth. element of Ri is given by 

RiM = Rsm^"' ( 4 ' 10) 

Then the beam power taking into account the effect of mean wind as prescribed by 
this model is 

P(/,u 0 ) = w + Riw (4-11) 

Equation 4.9 shows that the maximum magnitude of error introduced by this 
model would occur when the source was at an angle of 0° or 180° and the quantity 
(0 W - 0) was equal to tt. The direction of beam shift would depend on whether 0 was 
0 or 180 degrees. That is, the phase shift angle <pkn from Equation 4.9 is maximized 
for this set of conditions. The greatest change applied to the the velocity c occurs 
then since the cos0 factor is equal to il and the factor l/[c ■+■ V cos($w ^)] 
is maximized at l/(c - V ). Note that the cos0 factor scales down the effect as 0 
approaches ±tt/ 2 from either side. When 0 equals ±ir/ 2, 4> kn equals zero and there 
is no effect at all regardless of the magnitude of V because the cos# factor drives 
the whole term to zero. This is so because the plane waves are reaching all of the 
sensors at the same time (broadside). Also, there is no effect when u and p are 
orthogonal because in this case, the projection of the wind velocity in the direction 
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of u, cos(0w - 0), becomes zero. Figures 15 through 20 show representative 
beamformer outputs for frequencies of 30 , 45 and 90 Hz with 0 equal to 30 ° and 
60 °. The mean wind velocity V used was 13.4 m/s (44 ft/sec) and 0\v *as set equal 

to 6 for each case. 


C. Correction of Error Caused by Noise Model 1 

The effect of noise of the form of that proposed by Model 1 may be corrected. 
Since this effect is essentially a phase shift applied to the matrix Rs, (given in 
Equation 4 . 10 ), the matrix R x may be multiplied by another phase shift, synthesized 
from available wind data, that negates the phase shift caused by the wind. Note 

that =■ 1. 

Equation 4.9 defines the phase shift <t> kn applied to the fcnth element of the 
matrix Rs- If the respective elements of R x could be multiplied by a term e , 
then the effects of the mean wind would be completely eliminated. The exact values 
of <p kn are not known however because Equation 4.9 requires knowledge of the true 
bearing angle 0. This is the quantity that is sought. Instead, the assumed bearing 
angle d 0 is used in place of 8 to give the correction phase shift 7- 

(\ 1 \ ( 4 . 12 ) 

Tfcn = 27T/(n - fc)dcos0o^- - c + V cos(0w ~ 0o) / 

The fcnth element of the corrected spatial correlation matrix for the first noise model 
R\c is given by 

R\CM = R\,knS 


( 4 . 13 ) 
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EFFECT OF NOISE MODEL 1 



Figure 15. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 1 for / - 90 Hz an 

0 = 30°. 
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Figure 16. - Bartlett Estimate beamformer output comparing pure signal 1 <; a * e Wlt ^ 
output including atmospheric effects modeled by Noise Model 1 for / - 90 Hz and 


$ = 60°. 
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EFFECT OF NOISE MODEL 1 

1 .8288 m. U = 10. THETA = 30 DEGREES, f 
THETA_W = 30 DEGREES. V = 13.4 m/s 
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Figure 19. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 1 for / = 30 Hz and 
9 = 30°. 
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Figure 20. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 1 for / = 30 Hz an 
9 = 60°. 
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The beam power for the correction is then given by 

P(/,u o) = w + Ricw ( 4>14 ) 

At values of 0 O near 9, the correction is virtually exact. See Figures 21 through 
26 for representative results of these beamformer simulations for frequencies of 30, 
45 and 90 Hz and for values of 9 equal to 30° and 60°. The mean wind velocity 
was 13.4 m/s (44 ft/sec) and 9 w was set equal to 9 for each case. As these figures 
show, the beam shift caused by mean wind as represented by Model 1 can be totally 

eliminated. 


D. Noise Model 2 


The first proposed noise model assumed that the wind only affected the speed 
of propagation of the acoustic waves and not the direction. The second model to 
be shown now assumes that the wind only affects the direction of propagation and 
not the speed. As with the previous model, consider a wind at constant velocity V 
blowing over the array with bearing angle dw. Refer back to Figure 12. Let the 
wind vector and the true propagation vector add together to form a new propagation 
vector U 2 . The resultant i component is given by 


u 2x = — c cos 9 - V cos #w 


( 4 . 15 ) 


and the resultant y component is given by 


U 7 = — csin0 - V sinflw 


( 4 . 16 ) 
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CORRECTION OF NOISE MODEL 1 

c = 1.8288 m. M = 10. THETA = 3C DEGREES. I = 90 HZ 
THETA_W = 30 DEGREES. V = 13.* m/» 


7. CO 



Figure 21. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 1 and the case of the 
corrected Noise Model 1 for / = 90 Hz and 9 = 30 . 


CORRECTION OF NOISE MODEL 1 



Figure 22. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 1 and the case of the 
corrected Noise Model 1 for / = 90 Hz and 9 = 60°. 
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Figure 23. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 1 and the case of the 
corrected Noise Model 1 for / = 45 Hz and 0 = 30°. 
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Figure 24. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 1 and the case of the 
corrected Noise Model 1 for / = 45 Hz and 0 = 60°. 
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Figure 25, - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 1 and the case of the 
corrected Noise Model 1 for / = 30 Hz and 6 = 30°. 
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Figure 26. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 1 and the case of the 
corrected Noise Model 1 for / = 30 Hz and 9 = 60°. 
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A unit vector in this direction is u 2 with bearing angle d 2 given by 


9 2 = tan’ 




(4.17) 


u 2x J 

The propagation speed is also changed. The new propagation speed is c 2 and i 
given by 


= a / { 7 ?, "T U 2y 


(4.18) 


C2 = 

The maximum error caused by this model occurs when u and p are orthogonal 
regardless of the value of 9. 

This bearing angle is then used to calculate the adjusted spatial correlation 
matrix R 2 for the noise represented by the second model. The knth element of R 2 


IS 


(r*. — r n ) • U 2 , 

« 2 .fcn = S(f)S m (f) exp{j27r/ : } 


(4.19) 


Cl 


The beam power for this beamformer simulation is 


P(/,u 0 ) = W + R 2 w 


(4.20) 


Figures 27 through 32 give examples of output from the simulation of this noise 
model for frequencies of 30, 45 and 90 Hz and values of 6 of 30° and 60°. The mean 
wind velocity V was 13.4 m/s (44 ft/sec). In each case, 0w was set to 0 + 90° to 
obtain orthogonality of the wind and propagation vectors. 
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Figure 27. - Bartlett Estimate beamformer output comparing pure signal case wit 
output including atmospheric effects modeled by Noise Model 2 for / - 90 Hz and 

6 = 30°. 


EFFECT OF NOISE MODEL 2 



Figure 28. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 2 for / - 90 Hz an 

d = 60°. 
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Figure 29. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 2 for / — 45 z an 
9 = 30°. 
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Figure 30. - Bartlett Estimate beamformer output comparing pure signal ^case with 
output including atmospheric effects modeled by Noise Model 2 for / - zaD 
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Figure 31. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 2 for / = 30 Hz and 
9 = 30°. 


EFFECTS OF NOISE MODEL 2 
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Figure 32. - Bartlett Estimate beamformer output comparing pure signal case with 
output including atmospheric effects modeled by Noise Model 2 for / = 30 Hz and 
6 = 60°. 
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The first noise model was essentially a phase shift implemented by the multi- 
plication of a complex exponential. Thus the error it induced was easily corrected. 
The implementation of the second noise model is more mathematically complex. 
The error caused by noise of this form is located in the bearing angle 0 2 locating 
the unit vector -u 2 . In calculation of the beamformer, this angle is the argument of 
a cosine term. The cosine term itself is in the argument of a complex exponential. 
Since a linear array of equally spaced sensors is the subject of this analysis, the dot 
product in the exponential of Equation 4.19 can be written as 

(r fc - r n )-u 2 = {k - n)d cos 9 2 (4-21) 

9 2 is calculated from Equations 4.15 through 4.17. It is impossible to mathematically 
isolate 9 2 in order to correct for it. 

Efforts made to equate the effect of this representation of noise to an effective 
phase shift were unsuccessful. If this were possible, the same approach to correction 
as taken with the first model might yield favorable correction results. Another 
approach considered to correct for this error was the alteration of the steering 
vectors w t . Using the assumed bearing angle B 0 , it was possible to correct the 
beam shift at the true angle of propagation. However, this approach worsened 
the error in the “mirror image” of the beamformer output that exists due to the 
symmetry discussed in Chapter III. This symmetry must be acknowledged when 
using a linear array and the effects of such attempts at correcting for the error due 
to the wind would make the handling of the symmetry more difficult. Figure 33 
presents an attempt to use the steering vectors for correction. Note that at the true 
propagation bearing B = 45°, the corrected output and the pure signal output are 
coincident. However at the symmetric location of 9 0 = -45°, the uncorrected and 
corrected curves are coincident. The error is not corrected at this point. 
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Figure 33. - Illustration of an attempt to correct the Bartlett Estimate power 
output for the effects caused by Noise Model 2. Modification of the steering vectors 
Wi was used. The true bearing angle is as 45°. The symmetry of the linear array 
also causes an indication at -45°. At 0 O = 45°, the corrected plot coincides with the 
plot of pure signal output. However, at 6 0 = -45°, the corrected and uncorrecte 
plots coincide. This result complicates the symmetry problem because the ettect ot 
the correction is not symmetrical about 0$ = 0 • 
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CHAPTER V 

RESULTS, CONCLUSIONS AND 
RECOMMENDATIONS 

A. Results of Application of Noise Models to 
Bartlett Estimate Beamformer 

Figure 34 shows the behavior of the half power point beam width for the cases 
of pure signal, Noise Model 1, corrected Noise Model 1, and Noise Model 2. This 
figure shows that Noise Model 1 and Noise Model 2 yield a smaller beamwidth than 
the pure signal case and the corrected Model 1 case. For the set of parameters 
used in Figure 34, the data points for the beamwidth for the pure signal case and 
the corrected Noise Model 1 case are coincident. Also, the data points for Noise 
Model 1 and Noise Model 2 are coincident. The result shown is that both Noise 
Model 1 and Noise Model 2 have the same effect on the beam width for the same 
conditions at a particular frequency. 

Figures 35 and 36 show the error in true bearing angle (beam shift) for Noise 
Model 1, corrected Noise Model 1, and Noise Model 2. The error due to Noise 
Model 1 is generally greater than that caused by Noise Model 2. The correction 
for Noise Model 1 eliminates all of the error caused by Noise Model 1 but the cost 
of this correction is an increased half power point bandwidth as seen in Figure 34 . 
As all of the results show, neither of the noise models produce sidelobes. Also note 

that the beam shift is independent of frequency. 

Both noise models cause the beam to be more narrow and they both cause 
some beam shift. Since Noise Model 1 has the cos0 factor as in Equation 4.9, 
the error in true bearing location will depend on the value of 9 as well as the 
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Figure 34, - Relation between the half power point beamwidth and the frequency, 
/./ The data points for pure signal and Model 1 Corrected are coincident and the 
data points for Noise Model 1 and Noise Model 2 are coincident. The half power 
point bandwidth for the pure signal case and the corrected Model 1 case are always 
greater in this range of frequencies than for the cases of uncorrected Models 1 and 
2 at a given frequency. 
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Figure 35. - Results showing the error (beam shift) in the maximum of the Bartlett 
Estimate beamformer output for Noise Model 1, Noise Model 2 and the correction 
of Noise Model 1 versus frequency for 9 = 30°. 



Figure 36. - Results showing the error (beam shift) in the maximum of the Bartlett 
Estimate beamformer output for Noise Model 1, Noise Model 2 and the correction 
of Noise Model 1 versus frequency for 9 = 60°. 
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wind speed and direction. Noise Model 2 only depends on the wind speed and 
direction. Note in Figure 36 that, for 9 equal to 60°, Noise Model 2 causes more 
error. But, in Figure 35, Noise Model 1 causes the greater error in true bearing 
angle determination. 

Figures 37 and 38 show the effect of wind velocity V on the magnitude of 
the beam shift for Noise Models 1 and 2. In each case, the values of 9, d , M , 
0w and / are constant in each of the plots. These parameter values are not the 
same for both of the plots. In each of these plots, the wind bearing angle 0 W was 
selected to give the greatest amount of beam shift for the value of 9 used. It is 
clear that the error caused by each model increases with wind speed. The choice of 
22.4 m/'s (73.3 ft/sec) for V may seem high but it must be remembered that these 
are instantaneous velocities which may be sustained only for a fraction of a second. 
See Appendix B for discussion of atmospheric mean wind time scales. 

Figures 39 and 40 show the effect of the true bearing angle 9 on the magnitude 
of the beam shift over a range of 9. The values of M, d, f, V and 0\v are constant 
in each case (for each plot). Note that the maximum error for Noise Model 2 is seen 
in the middle of this range. However, for Noise Model 1, the maximum error occurs 
at the left end of the range for 9. This is so because the cos 9 factor in Equation 4.9 
goes to zero as 9 approaches tt/2. 

A maximum amount of beam shift for Noise Model 1 is shown in Figure 41. 
For a wind speed V of 13.4 m/s (44 ft/sec) and a value of 0\v selected to maximize 
4>kn from Equation 4.9, beam shifts of up to 10 degrees were obtained. This is 
observed as 9 is near 0 or 180 degrees and the quantity (0 W - #) is near tt. Note 
that as 9 becomes near 30°, the magnitude of the beam shift returns to the levels 
shown in the previous chapter - about four degrees and less. The wind conditions in 
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Figure 37. - Relation of beam shift to wind velocity for Noise Model 1 
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Figure 38, - Relation of beam shift to wind velocity for Noise Model 2. 
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Figure 40. - Relation of beam shift to true bearing angle 9 for Noise Model 2. 
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RANGE OF MAXIMUM BEAM SHIFT FOR NOISE MODEL 1 
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Figure 41. - Maximum beam shift observed for Noise Model 1 in the range of do 
near do equal to zero. 
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Figures 39 and 41 are not the same. As both figures show, the beam shift increases 
as 0 approaches zero degrees. It is not shown in Figure 39, but as Figure 41 shows, 
as 6 becomes less than about 20°, the amount of beam shift levels off and reaches 
a maximum magnitude of about 10 degrees. For wind speeds of up to 22.4 m/s 
(73.3 ft/sec), beam shifts near 0 and 180 degrees of up to 20 degrees may be seen. 


B. Conclusions Regarding Validity of 
Current Work 

Wind in the atmosphere degrades the performance of the Bartlett Estimate 
beamformer when used for tracking applications where the time for data sampling 
and processing is relatively short. This report has addressed one case of instanta- 
neous atmospheric conditions. This case is a wind of constant velocity blowing over 
an array of sensors in one direction. 
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Two models were proposed to mathematically represent the effect of mean wind 
on the Bartlett Estimate beamformer. Simulations were run and the effects were 
quantified. 

In the case of Noise Model 1, it was possible to correct, in closed form, for the 
effect of the mean wind. This model is the more idealized one. If in fact the mean 
wind may be represented by Model 1, the effect of the wind may be eliminated in 
the Bartlett Estimate algorithm. This was demonstrated in Chapter IV. 

For Noise Model 2, no closed form solution is available. The mathematics of 
the second model are more complex and an algebraic or trigonometric isolation of 
the error in bearing angle is not possible. More work is necessary to determine if 
perhaps an empirical correction for Noise Model 2 is possible. 

In the case of both of the noise models presented here however, the effects of 
the wind models are not as pronounced as expected. No obvious sidelobes were 
produced and no widening of the half power point beam width was caused. The 
only error caused by these models was a shift in the location of the maximum 
in the beam power which indicates the true bearing angle. The beam shift was 
generally four degrees and less. However, as was shown in Figure 41, beam shifts 
of 10 degrees are possible with Noise Model 1. It has been proposed that these 
beamforming techniques could be used to track targets 24.1 km (15 mi) away. A 
four degree error in bearing equates to approximately a 1.69 km (1.05 mi) error in 
location at a range of 24.1 km (15 mi). A beam shift of 10 degrees caused a 4.21 km 
(2.6 mi) error in location of the source. Therefore, this error could be appreciable 
depending on the conditions and the application. 

The work done in this report has three important results. 
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1) These models do cause one effect seen in actual results. This effect is the beam 
shift. Previously proposed models do not explain any of the three detrimental effects 
on the Bartlett Estimate beamformer. Thus, the present work is an improvement 
over prior work in this subject area in that at least one cause of beam shift has been 

identified. 

2) If the mean wind blowing over the array may be represented by the models 
derived and presented here, its effects have been quantified. And in the case of the 
first model, the effects of the mean wind may be corrected. If the amount of beam 
shift observed is deemed insignificant for the application, then effects of mean wind 
may be ignored and would be assumed to not effect the Bartlett Estimate in any 

crucial way. 

3) Since these models do not cause any sidelobes or increase in the half power point 
beamwidth as seen in actual experimental results, it must be concluded that the 
mean wind does not cause these effects. These effects must be caused by some other 
conditions in the atmosphere. There are many other instantaneous wind conditions 
which have not been addressed here. 

Some recommendations for continued work in this area are now made. 

C. Recommendations 

1) Conduct experimental verification studies in a controlled environment. The 
experiments would demonstrate the correctness of the two models presented here; 
experiments would also identify what atmospheric parameters need to be measured 


in addition to the sound. 
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2) Further literature search/experimental work is required to quantify the 
instantaneous wind speed and direction descriptors in the lower atmosphere. 
Address how to ascertain reasonable parameters for a wind gust and how to properly 
model turbulent eddies over a sensor array. 

3) Perform additional analytical studies (simulations) to determine the effects of 
array geometry, frequency, sensor spacing and the number of sensors deployed. 
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APPENDIX A 

DIGITAL SIGNAL PROCESSING CONSIDERATIONS 
The subject matter of this report is not concentrated on the practical imple- 
mentation of beamforming. However, two points should be mentioned. The relative 
magnitude of the beam power calculated is dependent on the number of FFT points, 
N, taken in the data sampling. As the length of the sampled input sequence is in- 
creased, more summations are carried out in the Fast Fourier Transform algorithm. 
This increases the magnitudes of the Fourier terms S(f ) in Equation 2.7. 

The choice of digital sampling frequency, f 3 , and N determine the frequency 
resolution of the beamformer. The frequency resolution is the “spacing” between 
frequency components and is given by 

A/ = — 

7 N 

If the sampled sensor output has frequency components that are closer together than 
A/, these components will not be represented as separate and distinct frequencies. 
Instead, “smearing” will result and the contribution to the spectrum from the 
in-between frequency components will be distributed among adjacent frequency 
components as discussed by Strum and Kirk (1988). The magnitudes of the Fourier 
Transforms of these components S{f) will be less than if the frequency component of 
the sampled data did not fall between Fourier components. When the beamformer 
calculates beam power for frequencies in a range where smearing has occurred, the 
magnitude of the beam will be somewhat less than those for frequency components 
which match exactly with frequencies in the sampled input sequence. Some digital 
signal processing techniques must be applied to maintain the performance of the 
beamformer. These techniques fall into the area of discrete spectrum analysis and 
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are beyond the scope of this report. A knowledge of the acoustic signature of the 
target is helpful. Selecting N and /, such that a fine frequency resolution is obtained 
over the range of interest aids in keeping a uniform beam power output. The 
sampling frequency,/, , must also be chosen such that frequency aliasing is avoided. 
f s must be at least twice the highest frequency in the sampled signal to avoid 

frequency aliasing. 

These issues were not addressed in the simulation done in this report. A value 
of unity for A / was used in all cases and all frequencies tested were integer values. 
These considerations are critical in practical application of the Discrete Fourier 


Transform to beamforming. 
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APPENDIX B 

ATMOSPHERIC MEAN WIND TIME SCALE 
Scales of time, length and velocity are commonly used in meteorology to 
quantify a time, length or velocity characteristic to a particular situation. A time 
scale is the length of time it generally takes for a certain action to occur as discussed 
by Stull (1988). There are many derived meteorological time scales. Some are for 
the stratification of the atmosphere. These time scales may be on the order of hours. 
Other time scales are for the viscous dissipation of energy by eddies on the order of 
one millimeter in size. Time scales such as these may be on the order of a fraction 
of a second. In this report, atmosheric phenomena are addressed which have time 
scales on the order of one second. The effect of mean wind on the performance of 
the Bartlett Estimate is the subject of this thesis. The motion of a large air mass 
is what causes the mean wind effects. According to Stull (1988), the time scale for 
such an action is x/U where x is the distance travelled and U is the mean velocity of 
the air mass. Depending on the relative magnitudes of these values, this time scale 
may be on the order of one second to one minute. For the cases shown in this report, 
x, the length of the array, is 16.5 m (54 ft) and L is 13.4 m./s (44 ft, sec). The time 
scale for this sensor arrangement is approximately 1.25 seconds. It is assumed 
that the wind conditions will remain constant over the array for 1.25 seconds. This 
is greater than the 0.5 second sampling time over which data is recorded. Therefore, 
the average of the data taken over the 0.5 second time window is assumed to be 
representative of the instantaneous atmospheric conditions at that time. 

For the case of a 22.4 m/s (73.3 ft/sec) wind as mentioned in Chapter V, this 
assumption still holds true. The mean wind time scale for the same array with a 
wind speed of 22.4 m/s (73.3 ft/sec) is about 0.75 second. While this length of 
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time is not an order of magnitude greater than the 0.5 second sampling time, it still 
considered to be of sufficient length over which to obtain representative average 

data. 

These assumptions are required for the Bartlett Estimate to be simulated 
as it was in this report. Atmospheric conditions are so unpredictable that such 
assumptions are necessary in almost any such study. The model presented in 
Chapter III that assumed the noise to be spatially white and homogeneous is an 
example of this necessity. In water, where that model is most often used, objects 
move much slower and longer sampling ang averaging times can be had. In the 
atmosphere, this is not the case. Therefore, these assumptions must be made to 
apply the technology which is already known to new and different conditions. 
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appendix c 


THREE MODELS FOR THE SPATIAL CORRELATION OF NOISE 
This appendix presents three models that were obtained from Burdic (1984) and 
essentially have been applied to three dimensional use in the ocean. The derivations 
of these models are not provided here and may be obtained in Burdic s text. See 

Figure 42 for illustrations of these models. 

1) Isotropic Noise - The isotropic noise model assumes that acoustic energy 
due to noise is distributed evenly over a sphere and uniformly over all frequencies 
and bearing angles. The ijth element of Q for isotropic noise is given by 


s in(47rdij/A) 

qij(da,\,T) = ( 4 V dij/X) 


(C.l) 


dij is the distance between the sensors, A is the wavelength corresponding to 
frequency /, and r is the time delay incurred over the distance traveled (r = d i} j c). 
These definitions also hold for the next two models presented. 

2) Semi-Isotropic Noise - The semi-isotropic noise model assumes that the 
acoustic energy due to noise is uniformly distributed over a symmetric portion 
of a sphere cut by an angle <t> 0 . The ijth element of Q for this model is given as 


qij(d{j , A, t) = sinc{2d^cos(0o)/A}cos(27rcr/A) (C.2) 


For an array on the ground, a value of zero degrees is used for <V The sine function 

is defined such that sinc(x) = [sin(x)]/x. 

3) Impulsive Noise - The noise is considered uniformly isotropic over a region 
everywhere in azimuth at an angle of elevation <f> 0 and zero elsewhere. The ijth 
element of the matrix Q is given by 


q t j(d t j,A) = cos(0 o ))Ml 


(C.3) 
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j 0 (.) denotes a first order Bessel function of the first kind. A value of zero degrees 

is used for © 0 for an array on the ground. 

For each of these three models, a simulation was run to examine the results 
when applied to the Bartlett Estimate beamformer of Equation 2.13. The simula- 
tions were done with a linear array of equally spaced sensors at a frequency of 90 
Hz with M = 10, d = 1.8288 m (6 ft), SNR = 1, and a true bearing angle 9 of 45°. 
The results of these simulations compared with beamformer output containing only 

pure signal are given in Figures 43 through 45 . 

In each case, there are no sidelobes present to indicate false sources. There 
are only extraneous levels of power output at assumed bearing angles other than 9. 
There is also no error in the indication of the correct bearing angle 9. Figure 43 
shows noise modeled as being isotropic as in Case 1 above. Here again the 
extraneous output is about 17 dB below the maximum value of beam power. 
Figure 44 shows the noise modeled as being semi-isotropic as described in Case 2 
above. This model yields the greatest level of extraneous power output at 11 dB 
below the peak value. Figure 45 shows noise modeled as impulsive as in Case 3 
above. There is a 13 dB difference between the main peak and the extraneous power 

i 

output caused by this model. The magnitude of the beamformer power output away 
from the true bearing angle 9 is not great enough to be considered detrimental to 
the beamformer performance. The absence of sidelobes indicate that these models 
do not cause the noise effects characteristic to actual beamforming results. Also, 
the half power point beam width of the main peak is not affected. Thus, these three 
models are inappropriate for tracking purposes. 
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Figure 45. - Bartlett Estimate beamformer power output for the case of pure 

signal compared with case including noise modeled as being impu sive a 

by Burdic (1984). 
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